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We have calculated the leading color-octet contributions to the production of J/ip 
particles in photon-proton collisions. Using the values for the color-octet matrix ele- 



ments extracted from fits to prompt J/ip data at the Tevatron, we demonstrate that 



distinctive color-octet signatures should be visible in J/ip photoproduction. How- 
ever, these predictions appear at variance with recent experimental data obtained 
at HERA, indicating that the phenomenological importance of the color-octet con- 
tributions is smaller than expected from theoretical considerations and suggested 
by the Tevatron fits. 
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The production of heavy quarkonium states in high energy collisions provides an im- 
portant tool to study the interplay between perturbative and non-perturbative QCD 
dynamics. A rigorous framework for treating quarkonium production and decays has 
recently been developed in Ref . [|1[] . The factorization approach is based on the use of non- 
relativistic QCD (NRQCD) to separate the short-distance parts from the long-distance 
matrix elements and explicitly takes into account the complete structure of the quarko- 
nium Fock space. This formalism implies that so-called color-octet processes, in which 
the heavy-quark antiquark pair is produced at short distances in a color-octet state and 
subsequently evolves non-perturbatively into a physical quarkonium, should contribute to 
the cross section. It has recently been argued in Refs. [^H that quarkonium production in 
hadronic collisions at the Tevatron can be accounted for by including color-octet processes 
and by adjusting the unknown long-distance color-octet matrix elements to fit the data. 

In order to establish the phenomenological significance of the color-octet mechanism it 
is necessary to identify color-octet contributions in different production processes. Color- 
octet production of J/ip particles in e + e~ annihilation and Z decays has been studied in 
Refs.0. In this note, we examine the color-octet contributions to J/ip photoproduction, 
7 + P — > J/ip + X, which proceeds predominantly through photon-gluon fusion. Elas- 
tic /diffract ive mechanisms || and reducible background processes can be eliminated by 
applying suitable cuts 0. According to the formalism developed in Ref.(l[], the inclusive 
production cross section can be expressed as a sum of terms, each of which factors into a 
short-distance coefficient and a long-distance matrix element: 

da( 7 + g -> J/tjj + X) = ]T d£( 7 + g^cc[n}+ X) (0|O J ^ [n]|0> . (1) 

n 

Here, da denotes the short- distance cross section for producing an on-shell cc-pair in a 
color, spin and angular-momentum state labelled by n. The NRQCD matrix elements 
(0\O J ^ [n]\0) give the probability for a cc-pair in the state n to form the J/ip particle. 
The relative importance of the various terms in (|l|) can be estimated by using NRQCD 
velocity scaling rules For v — *■ (v being the average velocity of the charm quark in 
the J/ip rest frame) each of the NRQCD matrix elements scales with a definite power of 
v and the general expression ([I]) can be organized into an expansion in powers of v 2 . 

At leading order in v 2 , eq.(|l]) reduces to the standard factorization formula of the 
color-singlet model ||. The short- distance cross section is given by the subprocess 

1 + 9^ cc[l, 3 S 1 ] +g (2) 

shown in Figfja, with cc in a color-singlet state (denoted by 1), zero relative velocity, and 
spin/angular-momentum quantum numbers 2S+1 Lj = 3 S\. Higher-order QCD corrections 
to the short- distance process @ were found to increase the color-singlet cross section by 
more than 50%, depending in detail on the choice of parameters P,|10 . 
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Figure 1: Generic diagrams for J/0 photoproduction: (a) leading color-singlet contribu- 
tion; (b) leading color-octet contributions; (c) color-octet contributions to inelastic J/ip 
production. 

Color-octet configurations are produced at leading order in a s through the 2 — > 1 
parton processes 



1 + 9 
1 + 9 



cc [8, ^o] 



cc [8, 3 P ( 



0,2J 



(3) 



shown in Fig.|T]b. The transition of the color-octet cc [8, 2S+1 Lj] pair into a physical J/ip 
state through the emission of non-perturbative gluons is described by the long-distance 
matrix elements (0\O J ^ [8, 2S+1 Lj]\Q) . They have to be obtained from lattice simula- 
tions or measured directly in some production process. According to the velocity scaling 
rules of NRQCD, the color-octet matrix elements should be suppressed by a factor of t> 4 
compared to the leading color-singlet matrix element ||11|| . Color-octet contributions to 
J/ip photoproduction can thus become important only if the corresponding short-distance 
cross sections are enhanced as compared to the color-singlet process. For the partonic 
cross section of the leading color-octet contributions we find 
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*(7 + g-+ cc [n] -> J#) = ?r<y(5 4m2 4mc2) Ei-^(7 + - cc [n])| 2 (O|0-W [n]\0) (4) 



with 



£|-M(7 + ^cc[8, 3 P ] 



(47r) 2 e 2 aa s 

4m c 
3(47r) 2 e 2 aa s 

4m 3 , 



El-M(7 + ^cc[8, 3 P 2 ])| 2 = M ^r S ■ (5) 

The partonic energy squared is denoted by s, and e c is the magnitude of the charm quark 
charge in units of e = ^Ana. We have checked that our results are consistent with those 
presented in Refs. [p|. |T2|1 . 

In the analyses Refs. [FJH , color-octet matrix elements have been fitted to prompt 
Jfip data from CDF [13| and found to be consistent with the NRQCD velocity scaling 
rules. The fit-values have, however, large theoretical errors due to the uncertainty in 
the heavy quark mass and the QCD coupling, unknown higher-order corrections, and 
higher twist effects [14j|, and should therefore mainly be regarded as order-of-magnitude 
estimates. For our numerical analysis we choose the color-octet matrix elements as listed 
in Table II consistent with the velocity scaling rules and the results obtained in Refs.|2]||. 
The color-singlet matrix element has been calculated from the J/ip wave function at the 
origin as obtained in the QCD-motivated potential model of Ref.||15[ and tabulated in 
Ref.||16||. At leading order in v 2 , the P-wave matrix elements are related by heavy-quark 
spin symmetry, (0\O J ^ [8, 3 Pj]|0) w (2J+ 1) (0\O J ^ [8, 3 P ]|0). 



(0|(9 J ^[n]|0> 


numerical value 


scaling 


(0\O J ^[l, 3 Si]\0) 
(O|0-W[8,%]|O> 
(0|O J ^[8, 3 P ]|0)/m 2 


1.16 GeV 3 
1(T 2 GeV 3 
1(T 2 GeV 3 


S 7 
3 7 



Table 1: Values of the NRQCD matrix elements used in the numerical analysis, with the 
velocity and mass scaling. 



Due to kinematical constraints, the leading color-octet terms will only contribute to the 
upper endpoint of the J/ip energy spectrum, z ~ 1 and p± ~ 0. The J/ip energy variable 
is defined by z — p ■ k^ / p • k 7 , with p, k^ a being the momenta of the proton and J/0, 7 
particles, respectively, and p± is the J/ip transverse momentum. The leading color-octet 



4 



3x10 



10° 



3x10 



10° 



3x10 



10' 



II 

- 7 + P J/V' 


111 

+ X 


1111 


1 1 1 1 1 1 
data: HI (prel.) — 


" o- [nb] 








_ z > 0.95, p ± < 1 


Ge V 




total 



















— — — ~~~ 


f8, 3 P ] 


- / 




.. — i 


— - 4 


^ [8,'So] □ 








1 






-1 


h 






^- 1 












1 








1 1 . 


-'"f ill 


1 i i i 1 


1 1 1 1 


1 i i i 1 i i i 



40 60 80 100 120 

Vs7 p [GeV] 



140 160 180 



Figure 2: Total cross section for J/ip photoproduction in the region z > 0.95 and 
p± < 1 GeV as a function of the photon-proton energy. Experimental data from [ |T7|j . 



and color-singlet contributions to the J/ip photoproduction cross section in the region z > 
0.95 and p± < 1 GeV are shown in Fig.|2|. The predictions are compared to experimental 
data obtained at HERA by the HI collaboration |TT|| . We have used m c = 1.48 ||15|| , 
a s (2m c ) = 0.3, and the values for the NRQCD matrix elements listed in Table |l| The 
parton cross section has been convoluted with the leading-order GRV parametrization of 
the gluon density in the proton [I8|, evaluated at a scale Q = 2m c . Relativistic corrections 
to the color-singlet channel fl9j enhance the large z region, but cannot change the order- 
of-magnitude suppression of the color-singlet process significantly. The large cross section 
predicted by the color-octet mechanism appears to be in conflict with the experimental 
data. This indicates that the value of the color-octet matrix elements (0\O J ^ [8, 1 S" ] 1 0) 
and (0|C J ^[8, 3 Pj]|0) is more than one order of magnitude smaller than expected from 
the velocity scaling rules and suggested by the Tevatron fits. It is, however, difficult 
to put strong upper limits for the octet terms from a measurement of the total cross 
section since the overall normalization of the theoretical prediction depends strongly on the 
choice for the charm quark mass and the QCD coupling. Moreover, diffractive production 
mechanisms which cannot be calculated within perturbative QCD might contaminate 
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the region z « 1 and make it difficult to extract precise information on the color-octet 
contributions. 



Diffractive processes can be eliminated by restricting the analysis to the inelastic 
domain z < 0.9 and p± > 1 GeV. Color-octet configurations which contribute to inelastic 
J/0 photoproduction are produced through the subprocesses 

7 + # -> cc[8, ^o] +g 
1 + 9 -> cc[8, 3 S!]+^ 

7 + <? -> cc[8, 3 P ,i,2]+^ (6) 

as shown in Fig{l]c. Light-quark initiated processes are strongly suppressed at HERA 
energies and can safely be neglected in a first step. The parton cross sections ([]) have been 
evaluated using the algebraic computer program FORM pg. Details of the calculation 



and analytic results will be presented in a forthcoming publication ||21|| . Here we merely 
note that the cross sections 7 + g — > cc [8, 1 5'o], [8, 3 -Po,2] + # are divergent for z — > 1 and 
— >■ 0, due to the g — > gg collinear splitting, Fig.[l]c2. This singularity can be absorbed 



into the renormalization of the parton densities via mass factorization ]2T|. In the inelastic 
region z < 0.9 and p± > 1 GeV the cross sections are well-behaved however. 

In Fig.|3| we show the color-octet and (leading-order) color-singlet contributions to in- 
elastic J/0 photoproduction compared to experimental data from HI |17|] obtained in the 



kinematical region z < 0.8 and p± > 1 GeV. For the central predictions we have used 
m c = 1.48 GeV and a s = 0.3, while the hatched error band indicates how much the 
color-singlet cross section is altered if m c and a s vary in the range 1.4 GeV < m c < 
1.55 GeV and 0.25 < a s < 0.35. Adopting the NRQCD matrix elements listed in Ta- 
ble g and (0\O J ^ [8, 3 5'i]|0> = 10" 2 GeV 3 we find that color-octet and color-singlet 
contributions to the inelastic cross section are predicted to be of comparable size. The 
short- distance factors of the [8, 1 S ] and [8, 3 -Po,2] channels are strongly enhanced as com- 
pared to the color-singlet term and partly compensate the O(10~ 2 ) suppression of the 
corresponding non-perturbative matrix elements. In contrast, the contributions from the 
[8, 3 <S'i] and [8, 3 Pi] states are suppressed by more than one order of magnitude. Since 
color-octet and color-singlet processes contribute at the same order in a s , the large size of 
the [8, 1 5'o] and [8, 3 Jb,2] cross sections could not have been anticipated from naive power 
counting and demonstrates the crucial dynamical role played by the bound state quantum 



numbers [p2j . Inclusion of the color-octet processes increases the inelastic cross section by 
about a factor of two, consistent with experimental data. This does, however, not proof 
the significance of the color-octet terms: as can be seen from Fig.[| the experimental data 
can be accounted for by the color-singlet process alone, once the theoretical uncertainties 
due to variation of the charm quark mass and the strong coupling are properly taken 
into account. This has also been demonstrated in the next-to-leading order analysis of 
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Figure 3: Total cross section for J/ip photoproduction in the region z < 0.8 and p± > 



1 GeV as a function of the photon-proton energy. Experimental data from |fPT |. 



Ref.||10||. No strong conclusions on the size of the color-octet terms can thus be deduced 
from the study of the total inelastic cross section. The same statement holds true for 
the transverse momentum spectrum, since, at small and moderate p±, both color-singlet 
and color-octet contributions are almost identical in shape. At large transverse momenta, 
P± 10 GeV, charm quark fragmentation dominates over the photon-gluon fusion pro- 
33LE3J. In contrast to what was found at the Tevatron [J25|, gluon fragmentation into 



cess 



color-octet states is suppressed over the whole range of p± [24 



A distinctive signal for color-octet processes should, however, be visible in the J/ip 
energy distribution da/dz shown in Fig.|4]. We have plotted color-singlet and color-octet 
contributions at a typical HERA energy of \/s 1P = 100 GeV in the restricted range 
p± > 1 GeV, compared to recent experimental data from HI ||26|| . Since the shape of the 
distribution is insensitive to higher-order QCD corrections or to the uncertainty induced 
by the choice for m c and a s , the analysis of the J/ip energy spectrum da/dz provides a clean 
test for the underlying production mechanism. The comparison with the experimental 
data, Fig.|, shows that the J/i/j energy spectrum is adequately accounted for by the color- 
singlet contribution. The shape predicted by the color-octet contributions is instead in 
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Figure 4: The J/ip energy distribution da/dz at the photon-proton centre of mass energy 
a/s 7P = 100 GeV integrated in the range p± > 1 GeV. Experimental data from [Z6fl. 



conflict with the experimental data. 

In conclusion, we have investigated color-octet contributions to the production of J/ip 
particles in photon-proton collisions. Color-octet processes would strongly affect the upper 
endpoint region of the J/ip energy spectrum and the shape of the J/ip energy distribution. 
However, these predictions appear at variance with recent experimental data obtained at 
HERA. We conclude that the values of the color-octet matrix elements (0|O J// ^ [8, 1 j Sq] |0) 
and (Q\O j ^ [8, 3 Pj] |0) seem to be more than one order of magnitude smaller than expected 
from the velocity scaling rules and suggested by the Tevatron fits. With higher statistics 
data it will be possible to extract more detailed information on the color-octet matrix 
elements, in particular from the analysis of the J/ip energy distribution in the inelastic 
region. 
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